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Abstract. We investigate the nature of a gravitational two-state system (G2S) in the 
simplest setup in Newtonian gravity. In a quantum description of matter a single motionless 
massive particle can in principle be in a superposition state of two spatially-separated 
locations. This superposition state in gravity, or gravitational cat state, would lead to 
fluctuations in the Newtonian force exerted on a nearby test particle. The central quantity 
of importance for this inquiry is the energy density correlation. This corresponds to the 
noise kernel in stochastic gravity theory, evaluated in the weak field nonrelativistic limit. 
In this limit quantum fluctuations of the stress energy tensor manifest as the fluctuations 
of the Newtonian force. We describe the properties of such a G2S system and present two 
ways of measuring the cat state for the Newtonian force, one by way of a classical probe, 
the other a quantum harmonic oscillator. Our findings include: (i) mass density fluctuations 
persist even in single particle systems, and they are of the same order of magnitude as the 
mean; (ii) a classical probe generically records a non-Markovian fluctuating force; (iii) a 
quantum probe interacting with the G2S system may undergo Rabi oscillations in a strong 
coupling regime. This simple prototypical gravitational quantum system could provide a 
robust testing ground to compare predictions from alternative quantum theories, since the 
results reported here are based on standard quantum mechanics and classical gravity. 


1. Introduction 

In this work we investigate the nature of a gravitational two-state system (G2S) in the 
simplest setup in Newtonian gravity. In a quantum description of matter a single motionless 
massive particle can in principle be in a superposition state of two spatially-separated 
locations. This superposition state in gravity or gravitational cat state would lead to 
fluctuations in the Newtonian force exerted on a nearby test particle. Though simple, this 
problem touches on three aspects, quantum, gravity and information. We give a brief sketch 
of the backdrop of quantum and gravity, then introduce the gravitational cat state and ask 
how to characterize and detect such a state in this setup. 
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Gravity and quantum are two cornerstones of modern physics which for a good part of the last 
century are the source of awe and wonderment in their amazingly wide range of applicability, 
and at the same time the source of persistent challenge in how to fuse them together into a 
unihed theory. Finding ways to unify quantum (Q) and gravitation (G), each having stood the 
test of time with great successes in their respective realms, is the goal of “quantum gravity” 
programs (schematically denoted as Q 0 G) , which has met with varying degrees of success 
(or dehciency, depending on one’s view and values). A more modest yet no less productive 
attempt is to place these two theories closely together (schematically denoted as Q © G), 
see what discrepancies this union may reveal, what difficulties we face in the interpretations 
and what new physical insights we may gain. Gonsidering that gravitation, as embodied in 
Einstein’s theory of general relativity (GR), is a classical theory for the large scale structure 
and dynamics of spacetime, and quantum held theory (QFT) for quantum matter down 
to the microscopic scale, this task is far from straightforward. The hnal analysis requires 
addressing the issues of quantum versus classical and micro versus macro which underlie 
almost all physical problems in contemporary physics. 

The theoretical effort of this vein can be traced back to the 60-70’s in quantum field 
theory in curved spacetime (QFTGST) [HEllS] which began in the late 60’s with cosmological 
particle creation studies and epitomized in Hawking’s 1974 discovery of black hole radiance 
1^. Focused efforts in seeking ways to regularize or renormalize the stress energy tensor of 
quantum helds made it possible to tackle the so-called ‘backreaction problem’ in hnding how 
quantum matter helds ahect the dynamics of spacetime. Solving the backreaction problem 
is at the core of semiclassical gravity theory (SGG) developed in the 80’s based on the 
semiclassical Einstein equation (SGE)[|] 

Following this vein, general relativists have probed the interplay between gravity and 
quantum largely from the angle of how quantum matter ahects spacetime (Q —?• G). Asking 
the question in the other direction (G—)■ Q), namely, how gravity could have an ehect 
on quantum phenomena, has been going on for just as long (e.g.. Ref. m) mainly by 
quantum foundation theorists. The foremost issue is why macroscopic objects are found 
sharply localized in space (their wave functions “collapsed” on dehnite locales) while those of 
microscopic objects extend over space. This contradiction is captured in the celebrated Gat 
of Schrodinger. One can very coarsely place these theories in three groups: The Ghirardi- 
Rimini-Weber - Pearle models [7j of continuous spontaneous localization, the Diosi-Penrose 
theories [HI [9] invoking gravitational decoherence, and the recent trace dynamics theory of 
Adler [10] which attempts to provide a sub-stratum theory from which quantum mechanics 
emerges. A nice description of these theories can be found in a recent review by Bassi et al 

[m. 

These efforts spaning several decades are now called back to attention mainly due to the 

f For references in these developments, see, e.g., the bibliographies in the monographs listed above and 
references listed in a recent overview [5] with relevance to the issues addressed in this paper. 
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fast development of experimental capabilities to put these theories to tests. In fact, a new 
name “gravitational quantum physics” (see, e.g.. Ref. [12] ) is coined by experimentalists to 
capture the convergence of interest from both the quantum foundation and the quantum + 
gravity communities. The interesting array of ongoing experiments using atomic-optical, 
molecular spectroscopy, nanomechanical and other schemes to test alternative quantum 
theories and/or the effects of gravity on quantum systems also bring in lively issues in the 
quantum-classical and micro-macro interfaces. In this context, issues related to quantum 
information such as quantum decoherence and entanglement arise. This last aspect is the 
target scope of our present paper. 


1.2. Quantum Superpositions involving Gravity 


We now add quantum information considerations to quantum systems involving gravitational 
interaction. A simple problem of foundational interest is the superposition of two 
macroscopically distinct states. This simple yet basic state, which goes by the name of a “cat 
state”, has been studied and tested out for some nongravitational systems [I3] . Experimental 
schemes for testing quantum superposition of mechanical objects (e.g., mirrors, cantilevers) 
by their interaction with light and electric devices have also been explored na El uni nTj in 
the frontline research areas of opto- and nano-electro-mechanics. 

As illustration of a gravitational cat state, consider the quantum description of a 
stationary point mass m localized around x = 0 with spread a, described by a Gaussian 
wave function with zero mean momentum. 

The position x of the particle is a random variable described by the probability distribution 
|'0o(x)p. However, according to Newton’s law, a probability distribution for x, dehnes a 
probability distribution for the Newtonian force acted on a particle of mass mo located at R 

Gmmo 


F = 


IR 


R 


( 2 ) 


For |R| » a the fluctuations of the Newtonian force are negligible, which leads one to view 
it as a deterministic variable. 

Now consider a cat state, i.e., a superposition of two Gaussians, each located at ±^L 
and with zero mean momentum. 


V’(x) 


1 1 
71 (27ra2)3/4 


_(x-L/2)2 _(x-L/2)2 

e -|- g 


(3) 


If L is of the order of magnitude of R, the fluctuations of the Newtonian force ([^ are non- 
negligible. Since the force is a function of x, and x is described by an operator in quantum 
mechanics, the Newtonian force should also be described as an operator. But then, so would 
be the gravitational potential. In this sense, the cat state for the point mass has generated 
a cat state for the gravitational held. 
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What is the right theory to use for the description of a gravitational cat state? What 
characteristics do they possess and how they can be probed are the questions we ask in this 
paper. These inquires are the starting point of any investigation into quantum informational 
issues of gravitational systems. 

1.3. What is the right theory to use? 

In answer to the first question above, we may begin searching for a theory which describes 
quantum matter existing / propagating in a classical gravitational held, or QFT in curved 
spacetime. Further demanding self-consistency in the solutions of the dynamical equations 
for both classical gravity and for quantum matter we reach the semiclassical gravity theory 
described by the semiclassical Einstein equation (SCE), (also known as the Moller-Rosenfeld 
equation [18] in the alternative quantum theories community): 

G,, = SttG < T^,(<F) > (4) 

where is the stress energy tensor of the matter held, here represented by a scalar held 

<h, and <> denotes taking the expectation value with respect to a certain quantum state, 
usually the vacuum state. 

However, we run into a problem immediately. Consider a single particle. For a 
superposition of two states each localized at x = ±|L, the SCE eqn makes the prediction that 
the particle is most likely to be found at the center x = 0, rather than at either side. This 
intrinsic dehciency of SCG has been pointed out in e.g., HHEO]. Some authors [2T1122| use 
this as an indication that gravity needs be quantizec[|j One needs to include considerations 
of the quantum correlations, or equivalently, quantum huctuations of the mass density of the 
system to be able to capture the physics of the superposition state. 

There is such a theory, known as stochastic gravity introduced in the 90’s whose main 
challenge was to hnd a lawful place for the fluctuations of quantum helds as a part of the 
total source driving the Einstein equation. Adding the expectation value of the stress energy 
bitensor, also known as the noise kernel, to the expectation value of the stress energy tensor 
upgrades the SCEq to the Einstein-Langevin eguation m, which enables one to solve for 
the induced metric fluctuations (the poetic ‘spacetime foam’) in a self consistent manner. 
Stochastic gravity [33] has been developed and applied to strong field situations such as 
structure formation in the early universe [35] and black hole fluctuation and backreaction 
issues [51]. Here, in the opposite regime, of weak held and nonrelativistic limit, stochastic 
gravity based on GR -|- QFT is also the theory suitable for addressing quantum information 
issues at laboratory accessible low energy ranges [j]] 

§ Note the subtle yet important distinction between quantizing the weak perturbations of gravity (gravitons 
as spin 2 fields) and quantum gravity, defined as theories for the microscopic structures of spacetime. The 
former is on the same footing as phonons in regard to the latter as atoms [53] . The former can exist at 
today’s low energy while the latter pertains to Planck scale physics which we have limited facts to speculate 
on. 

II In terms of theoretical structure, stochastic gravity occupies the intermediate range as one progresses 
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For readers familiar with the essentials of stochastic gravity theory, the main character 
in our exposition, the mass density correlator, is just the weak field nonrelativistic limit of the 
noise kernel. However, our presentation below does not assume any familiarity of stochastic 
gravity. It should be understandable with some minimal knowledge of nonrelativistic 
quantum helds and classical gravity. 

1.4- Gravitational Cat State 

In non-relativistic gravity, a mass density distribution /i generates a gravitational potential 
0 via the Poisson equation 

VV = -dTrGp. (5) 

A fluctuating mass density implies a fluctuating potential and thus a fluctuating force on 
any test particle. Hence, a quantum state that corresponds to a superposition of states 
with different mass density distributions gives rise to a ‘superposition’ state of gravitational 
potentials, i.e., a gravitational cat state. 

Gravitational cat states are not new. The reason they are not in our conscious state 
of mind, much less the object of serious study, is likely because they have been viewed as 
undesirable, as targets of elimination, from the start. This is the spirit of many alternative 
quantum theories Ini- 

Gravitational cat states are rarely mentioned in their proper setting, from quantum 
field theory in curved spacetime to quantum gravity, because these theories mostly aim at 
describing ultra-high energy processes near the Planck scale. Only in the field of quantum 
cosmology in the context of the “wave function of the universe” EH or the “universal wave 
function” (many-worlds interpretation) [38] is an analogous problem discussed, because the 
generic solution to the Wheeler-DeWitt equation is a cat state, i.e., a superposition of two 
distinct cosmological spacetimes—see, for example. Ref. [39]. Our aim here is to explore 
the nature and consequences of gravitational cat states in the low energy realm that is more 
amenable to laboratory experiments. 

Gravitational two state (G2S) systems should play a prominent role in the new held called 
“gravitational quantum physics” [12], ushered in by leading experimentalists working on 
atomic-molecular, opto- and electro-mechanical schemes wanting to explore issues involving 
quantum, information and gravitation. 

1.5. Gravitational Decoherence 

It is perhaps useful to recall Penrose’s inhuential argument against the existence of 
gravitational cat states [9], that is employed in support for gravitational decoherence. It 

(‘bottom-up’) from semiclassical to quantum gravity, in the sense that when higher correlation functions of 
quantum matter helds are added as source of the Einstein equation the quantum coherence of the gravity 
sector begins to appear and increase in degree with the order of the correlation hierarchy. See description in 
|36j . The SCG as a mean held approximation capturing only the lowest order is understandably inadequate 
in this regard. 
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proceeds as follows. 

In quantum theory, time is an external parameter manifested in the evolution of the 
quantum states through Schrodinger’s equation. In general relativity, time is part of the 
spacetime structure, which is dynamical and depends on the conhguration of the matter 
degrees of freedom. This contradiction is manifested clearly when considering superposition 
of macroscopically distinct states. Each component of the superposition generates a different 
spacetime. Since there is no canonical way of relating time parameters in different spacetime 
manifolds, there is a fundamental ambiguity in the choice of the time parameter in the 
evolution of the quantum state. This ambiguity is manifested even at low energies, when the 
gravitational interaction can be effectively described by the Newtonian theory. It suggests 
a physical mechanism of gravity-induced decoherence for superpositions of macroscopically 
distinct states. 

We agree that the ambiguity in the dehnition of the time coordinate in a cat state is 
real, as far as our present knowledge of fundamental physics is concerned. In fact, we have 
rephrased this argument in terms of gauge independence for the symmetry of space and time 
reparameterizations in Ref. BD]. However, the conclusion that a cat state must decohere 
does not follow. As we have argued in Ref. [H], strong additional assumptions about the 
character of such ambiguities or fluctuations are needed in order to infer that they are agents 
of decoherence. Schrodinger’s cat may very well coexist peacefully with temporal ambiguities 
or fluctuations. 

1.6. The present study 

We study the stress energy tensor and its fluctuations at the level of a single non-relativistic 
particle. In the non-relativistic limit, the dominant contribution to the stress-energy tensor 
comes from the mass density /i(r). The corresponding quantum operator is proportional to 
the number density operator that dehned in terms of the particle creation and annihilation 
operators, or equivalently, the non-relativistic quantum helds. Thus, expectation values and 
correlations of the mass density can be dehned as expectation values of products of /i(r,f). 

1.6.1. Properties of fluctuations. With the above characterizations, we establish two 
important properties of the mass density huctuations. 

First, mass density fluctuations are very strong in single-particle systems. They are of 
the same order of magnitude as the mean mass density. This effect, though somewhat unusal, 
is in line with the earlier Endings of Ford, Hu and their co-workers [251 [261127|, where the 
huctuations of the energy density of quantum helds in Minkowski space, Casimir geometry, 
Einstein universe are found to be of the order of their mean values. This was one of the 
reasons why huctuations need be included in the consideration of the inhuence of quantum 
matters on spacetime structure and dynamics, as highlighted in stochastic gravity theory. 

This result is in contrast to the nature of mass-density as dehned in classical 
hydrodynamics, where huctuations are usually assumed negligible. However, the suppression 
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of fluctuations in hydrodynamics is statistical. It arises out of the large number of particles 
constituting the physical system and quantum fluctuations are usually unimportant in 
comparison to statistical fluctuations. However, when considering a single- or few-particle 
system statistical fluctuations on the number of particles in a given region, while important, 
do not, in general, dominate over quantum meachanical fluctuations. 

Second, there exists a relation between the correlation functions of the mass density 
and position histories of the associated quantum particle. The correlation functions of the 
mass density correspond to values of the decoherence functional for position histories of the 
particle. The decoherence functional is a complex functional of pairs of histories that appears 
in the consistent/decoherent histories approach to quantum mechanics [2H112SI EHl EH E2]- 

1.6.2. Probing a Gravitational Two-State (G2S) System. We consider a particle prepared 
in a cat state associated with the two minima of a conhning potential and we calculate 
the corresponding gravitational force. Two cases are considered. The hrst corresponds 
to a classical probe that records continuously the gravitational force. We compute the 
full probability distribution for the force measurements, and we show that they satisfy an 
exponential law and are distinctly non-Markovian. The second corresponds to a quantum 
probe, an harmonic oscillator acted upon by the Newtonian force. We show that the system 
is described by the Jaynes-Cummings Hamiltonian, and we determine the probe’s evolution 
in the relevant regime. We find that tunneling between the two minima of the potential leads 
to Rabi oscillations for the probe. 

1. 7. Organization of this paper 

We will hrst describe in Sec. 2 the stress energy correlations in nonrelativistic held theory, 
and we show their relation to the decoherence functional and to the Wigner function that 
describe the time evolution of a single particle. In Sec. 3 we consider a quantum system 
characterized by a conhning potential with two local minima such as a double well potential. 
In Sec. 4 we present two ways to probe the properties of such a gravitational two level 
system. In Sec. 5 we summarize the important features brought out in our analysis of this 
prototypical gravitational quantum system. 

2. Stress-energy correlations 
2.1. Nonrelativistic N Particle System 

Consider a scalar quantum held 0(r) and its conjugate momentum 7r(r) expressed in terms 
of the creation and annihilation operators and al 



( 6 ) 


( 7 ) 
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For a free field, the Hamiltonian operator is 


H = 


(Pk 


t. 


(27r 




( 8 ) 


where Wk = + 




In the non-relativistic regime we dehne the helds 


^(r) = 


d^k 


j (27r)3 
and to leading order in |k|/m, 


Cl]iC 


ikr 


p\r) = 


d^k 


t —ik-r 


(2vr) 




0(r) = 


^/2m 

The Hamiltonian then becomes 


■^(r) + '^’*'(1 


7r(r) = — 


■^(r) — '^’*'(1 


H = m / dr-^^(r)-^(r) — -— / dvip\ir)V‘^p{v). 


(9) 


( 10 ) 


( 11 ) 


We will denote the second term in Eq. (11) as Ho because it corresponds to the Hamiltonian for 
N non-relativistic particles. The hrst term in Eq.([TT|) corresponds to iVm, for an iV-particle 
state. Hence, the number operator N is 


N= dvp^ {ir)'ip{v) 


( 12 ) 


This suggests that m'^i(r)'^(r) can be identihed as the mass-density operator 

We include the effect of a conhning potential V (r) , by modifying the held Hamiltonian 


H = m / dr'^^(r)'^(r) -f / dvp^(v) 






(13) 


2.2. Mass-density correlations, noise kernel 

In the non relativistic limit, the dominant component of the stress tensor T^i, is the energy 


density, which is dominated by the mass density, namely 

T^^(r,f) = 5j5Xr,f) (14) 

Thus, it suffices to calculate the correlation functions of the Heisenberg-picture operator 

/i(r,f) = e'^'/i(r)e-'^h (15) 

We assume an one-particle state 

10) = y dr0(r)^/>t(r)|O), (16) 

where 0(r) is the one-particle wave-function. 

We hnd 

(/i(r,f)) = (0|/i(r,f)|0) = m0*(r,f)0(r,t) (17) 

(/i(r,t)/i(r',T)) = 0*(r,f)0(r',F)G'(r,t;r',T) (18) 
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where 4>{r,t) is the time-evolved single particle wave function and G(r,t;r',t') is the one- 
particle propagator, 

= (r'|e-*^(*'-*)|r). (19) 

For a free particle, 


m \3/2 




?m(r —r'l , , 

We note that the two-point correlation function is complex valued and cannot dehne a 
stochastic process. However, the real part, 

^(r, t; r', F) = Re{fi{r, t)p(r', t')), (21) 

known as the noise kernel, corresponds in some cases to the two-point correlation function 
of a stochastic process. 

Of importance is also the connected two-point correlation function for the mass densities 


.02 


i;(r,«; r',«') = {/i(r, t*)) - (/i(r, i)){;i(r', i')). 


( 22 ) 


2.3. Relation to the decoherence functional 

Let us ponder for a moment the meaning of the mathematical expressions above. Note hrst 


that the expectation values (18) do not correspond to the correlation functions of a physical 


process, because in realistic systems the mass density is not dehned at a sharp spacetime 
point but smeared over a hnite spacetime region. In actual experiments, the particles under 
consideration (atoms) have a hnite size d and it is meaningless to talk about mass densities 
at scales smaller than d, unless one has a detailed knowledge of the particle’s internal state. 

For this reason, rather than the exact mass density function, we consider a smeared 
mass density function (for more about applying smearing function for the evaluation of 
distributions, see, e.g., |271E1]) 


fis{r,t) = j dr'f{r-r')fi{r',t), 

for some smearing function /(r) of dimension [length]"^, 
smearing scale £ is dehned by the condition £^ = l//(0). 

We dehne the positive operator 

Pr= f - r')|r')(r'|, 


(23) 


centered around r = 0. The 


(24) 


where ^'(r) ;= /(r)//(0). The operator P^ represents a sampling of position with a width Sx 
around r. 

If the sampling function is a characteristic function of some set (i.e., if g‘^ = g), then P^ 

_ r2 

is a projection operator. Here, we will consider Gaussian functions of the form g{r) = e ^ 
in which case Pr is an approximate projector. The corresponding smearing function is 

1 

-^— rr 


/(r) = 


(27r)3/2c 


(25) 
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and corresponds to i = \/^Sx- 

The correlation functions of the mass density become 

HI 

{Ps{r,t)) = -{<P\PxM (26) 

9 

7Tl '' 

{/i,(r,t)/i,(r',t')) = -^{<P\PrtPr't'\(i)), (27) 

where Prt = is the Heisenberg-picture evolution of Pr- 

The expectation value of the smeared mass density is proportional to the probability of 
a position measurement at time t. The two point correlation function is proportional to the 
decoherence functional 

T>(r,t;r',T) := (0|Prf-Pr't'|0) (28) 

for a pair of histories one corresponding to a position record r at time t and the other to a 
position record r' at time t'. As explained by the decoherent histories approach to quantum 
mechanics, the vanishing of the off-diagonal elements of the decoherence functional in an 
exhaustive and exclusive set of histories is a necessary condition for assigning probabilities 
to this set of histories. The assignment of probabilities implies that the time evolution of 
the system can be expressed in terms of classical equations of motion, possibly including a 
stochastic component. 

Typically an approximate projector remains close to a true projector under time 
evolution, so P^^ ~ Prt- Of course, the equality is exact if Pr is a sharp projector. Then, 
(/rs(r, t)/^s(i', t)) = p(/rs(i', t)), and the connected correlation function becomes 

7Tl 

r/(r,t;r,t) = —(/r,(r,t)) - (/i,(r,t))2 > 0. (29) 

For t = t', P(r, t; r', t) ~ 5(r — r'); thus, the correlation function (/is(r, t)/is(r', t)) vanishes at 
equal times, unless r = r' within an accuracy of order i. 

The values of the decoherence functional P(r,t;r',f') for general arguments can differ 
signihcantly from zero only if the vectors Prt\4>) and Pr't'\4>) overlap substantially. Consider 
as an example the case of a free particle and assume that |0) is well localized in position 
around some point fq. The vector Prt\(p) is supported only on the momentum components 
p of 10) such that p = m(r — ro)/t, and similarly Pr't'|0) is supported only on the 
momentum components p of |0) such that p = m{r' — ro)/t'. The two vectors overlap if 
(r — ro)/t — (r' — fq) /t'. If this condition is satished, the values of the decoherence functional 
are determined by the time-of-flight momentum 

f' - f' 

Ptf = ^ ■ (30) 

as 


P(F,f;F',F) ~ |(0|prF)P^[(r - Fo)/t - (f' - Fo)/t']. (31) 

The associated connected correlation function is very different from a delta function on t and 
f; hence, fluctuations exhibit highly non-Markovian behavior. 

We note that Eq. (31) is relevant for well localized initial states, and does not apply, for 
instance, to a cat state. 
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Next, we proceed to evaluate the smeared correlation functions (27), using the Wigner-Weyl 


representation, in which an operator A on the Hilbert space of a particle is represented by a 
function on the associated state space. 


y \ Ap y 


^i(x,p) = y dy(x--|H|x+-)e 

We denote the Wigner-Weyl transform of Pr,i as Pr,t- For a free particle, Pr,t is 

Fr,t(x,p) =^(r-x - —). 

m 


(32) 


(33) 


The Wigner-Weyl transform of the product PrtPr't' that appears in Eq. (27) is 

approximated by 


Fr,ty.f{^,P) ^ ■f’r,t(x,p)fr',t'(x,P). 


(34) 


The approximation (34) is meaningful in a quasi-classical regime where x and p correspond 


to coarse-grained observables, like, for example, the position and momentum of the center of 
mass of a composite particle. 


With a smeared mass density function (|23|) this yields, 

1 / r -h r 


Fr,t;r',t'(x,p) = OXp 

(f - t'f 


AnPsI 


p — m 


s 

r — r 
t-t' 


X — 


X 

/\ 2 ' 


P(t + t') 
2m 


Substituting into Eqs. (27), we obtain 


(ft(r,i)) = ^ j ^^'ro(x,p)Fr,,(x,p) 

{/j.(r,i)/j,(r',i')> = ^ j ^^'ro(x,p)f;j;...j.(x,p) 


(35) 

(36) 

(37) 


where Wq is the Wigner function of the initial state. 

For scales of observation much larger than s^, we substitute the Gaussians in Eq. (33) 
by delta functions 


Fr,t(x,p) ~£^5(r-x-^' 
Fr,t;r',t'(x, P) — ^^<5 ( X — 


m 
r -I- r' 


2m ^ m^ ” 


whence, with the use of (35), 
(hs(rG)) 


f dpWo{r - —,p) 
(27rj-^ J m 

m^ r- _l_ r' 

Wn( —T- 


(27r)3(f — t'y 


(r-r')(t + t') r-r' 

, m- 


2{t - t') 


t-t' 


(38) 

(39) 

(40) 

(41) 
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The two-point correlation function is, therefore, straightforwardly calculated for any initial 

(r-i-O(w') j ^ m— 
2(t-t') P ^ t-t' • 


state by computing its Wigner function and setting x = 


Special cases Next, we consider an initial wave-function 'i/jq with vanishing mean momentum. 
We assume that the momentum spread is so small, that at a macroscopic level of observation, 
the p-dependence of the Wigner function is essentially that of a delta-function, lTo(r,p) = 
(27r)^5(p)|'0o(i’)P- Then, we obtain 

(/i,(r,t)) = m|'0o(r)r- (42) 

{/i^(r,t)/i^(r',t')) = mVo (r) - r). (43) 

Both the expectation value of and the two-point correlation function of fig are time- 
independent for this class of initial states that includes both the Gaussian state Q and 
the cat state (1^. Both mimic stationary states within this approximation, because of the 


vanishing mean momentum. We also note that the delta function in Eq. (43) is meaningful 


for variation of the arguments at scales much larger than £. For r ~ r', it is appropriate to 


substitute the delta function by the smeared delta function /(r) as in Eq. (25) 


The expectation value t)) is equivalent to a classical mass density p(r) = m\'ipo{r)\'^. 
However, the two-point correlation function does not correspond to fluctuations of the 
classical mass-density p(r), because the latter would be proportional to |V’o (i") 1"^ rather than 
to (r) P- Thus, the correlation function (43) reveals a characteristically quantum property 


of the density fluctuations, with no classical analogue. 

We estimate the relative size of the fluctuations by evaluating the ratio 

virG;rG) ^ I . , 

^ ^ (/r.(r,f))(/r.(r,t)) ^ ^ 

C{r) is not restricted to values much smaller than unity, thus, the correlations are of the 
same order of magnitude of the mean value, or even larger. 


2.5. Correlation functions for a stochastic process 

The quantum correlation functions of the mass-density are, in general, complex-valued, 
and they do not correspond to the correlation functions of a stochastic process. Such a 
correspondence would be highly desirable, because when considering the measurement of 
such correlations, it is much more convenient to treat them as subject to the rules of classical 
probability theory 

In general, the quantum correlation functions dehne a stochastic process only if specihc 
decoherence conditions are satished In this case, it can be proven that, for Gaussian 
processes, the real part of the quantum two-point correlation function (the noise kernel) is 
the two-point correlation function of a stochastic process. This implies that the observable is 
essentially classical and its evolution can be modeled using classical probability theory—see 


Refs. [32| I43| 142] . This is the case of the two-point function (43), albeit in a rather trivial 
sense, because the correlation function is time-independent. 
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The correspondence of the correlation functions to the decoherence functional suggests 
of a different way to dehne statistical correlation functions associated with the mass density. 
We dehne 


( TTl \ ^ 

—j P„(ri,fi;r2,f2; • ■ ■ 


where 


Pni^ly ti, r 2 , ^ 2 ) • • • J ^ntn) Tv Pr^tn ■ ■ ■ .^-2*2.^-2*2 ■ ■ ■ Pi 


T* 'f' 

i n f'n 


(45) 


(46) 


is the probability that n position measurements at times ti, i = 1,2,... ,n yield values rj. 


The time-instants are ordered: ti < t 2 < ... < tn- Eq. (46) applies if the operators Prt are 


sharp projectors. If P^t are approximate projectors, we substitute each P^-ti in Eq. (46) with 


P 


Viti ■ 


The probabilities (46) are well dehned and they may be employed in order to calculate 


the response of macroscopic detectors to the quantum system. However, they do not dehne 
a stochastic process, because they do not satisfy the Kolmogorov additivity condition that 
is essential for the dehnition of a stochastic process. 


Pn-iii 2 yt 2 ;---',intn) ^ (r,, f,; r 2 , ^ 2 ; ■ ■ ■ i r„i„). 


(47) 


Furthermore, the probabilities (46) are strongly dependent on the coarse-graining scale i 


In general, £ is a scale that is determined not only by the characteristics of the 
measured system (e.g. the size of an atom) but also from the method of sampling and 


precision of measurement of its position. This implies that the probabilities (46), and. 


consequently, the statistical correlation functions (45) depend on the specihc set-up by which 


the mass density is being measured. In particular, they depend both on the sampling length 
i, but also on ‘how often’ the system is being sampled, i.e., the number n of measurements 


as in (46) in a given time-interval. Hence, unless we consider systems in which the mass 


density behaves classically (as in (43)), different measurement schemes will lead to different 


predictions for the fluctuations of the mass density, and, consequently, of the associated 
Newtonian force. 


3. A gravitational two-state system 

In this section, we consider a quantum system characterized by a conhning potential with 
two local minima located at r = ±^L. We label the minima as + and —. At the macroscopic 
level, we can only distinguish whether the particle lies in the -|- region or in the —region. We 
take the direction of L to correspond to the x axis. Then, the projection operators onto the 
± regions are given respectively by 


P+ = 


P+ = 


dx 


dx 


’ OO 

poo 


dy / dz\x,y,z){x,y,z\ 

J OO 
poo 

dp / dz\x,p,z}{x,p,z\. 


(48) 
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By this definition, + P_ = 1. We will employ indices a, b, c and so on, that take two values 
± in order to denote the two minima. 

We consider an initial state 


IV') = C+I+)+c_|-). (49) 

where |+) and |—) are states localized around the minima + and — respectively, i.e., they 
satisfy 

Pa\a) = Sab\a), a,b = ±. (50) 


Thus, the system can be expressed as a quantum two state system (or qubit), where 


i+) = 




(51) 


In this approximation, it makes no difference whether we consider the projectors P± or any 
other projector that satisfies Eq. (50). The smearing scale i is identified as the width of the 
finest projectors that satisfy Eq. (50). Hence, it is of the order of magnitude of the width of 
the localization region around the potential minima. 


3.1. Trivial dynamics 

First, we assume that the tunneling probability between the two minima is so low, as to be 


negligible. This implies that UlPaUt = Pa- 


j - < , 1 1 2 

{Ts{a,t)) = \Ca\ 

(52) 

2 

TTl 

{frs{a, t)Ps{b, t }) 

(53) 


for all a,b = ±. Thus, the correlation functions are time-independent. A measurement at 
any moment of time (through the exerted gravitational force) will find the system in one of 
the minima with probability given by |cap. The Sab in the two-point correlation function 
signifies that the measurement collapses the wave-function completely and any subsequent 
measurement will find the system in the original minimum. 


3.2. Tunneling dynamics 


We consider a two-level system, in which tunneling effect is an important contribution to 
its dynamics. A double-well system is typically characterized by a time-scale e at which the 
transition amplitudes stabilize. Hence, we can describe the system’s evolution in terms of 
successive time-steps. For a two-state system |±) the evolution operator u for each time-step 
e is 


u = 


cos 6 sin 
— sin 6e~^^ cos 6 


( 54 ) 


for some angles 6 > 0 and y, where 6^ << 1. 
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For t » e, we take the continuous-time limit by computing vP and setting n 
obtain an one parameter family of unitaries 


Ut = 


cosf 

sinfe-*>^ 


sin 
cos y 


for V = 20/e, which correspond to a Hamiltonian 
H = //(cosydi -1- sinxd 2 ). 

Then, the evolved projectors P±t = U^P±Ut are 


P±t ~ 2 ^ 


St = 


cos ut 
sin 


sin 

— cos ut 


tie. We 


(55) 


(56) 

(57) 

(58) 


where a = ±1. 

The quantum correlation functions are 

7Tl 

{/i(a, t)) = —— [1 -|- a{S cos lot + (3 sin vt)] ( 59 ) 

2 

Tfl 

{fi{a2, t2)fi{ai, ti)) = — [1 + ( 5 (ai cos uti + 02 cos ut2 + 0102 cos u{t2 - ti)) 

+I 3 {ai sin vti + 02 sin iJt2) — *70102 sin u{t2 — ti)]. ( 60 ) 

where 0,01,02 = ±. The real numbers 5,7 and (3 describe the initial state: 5 = |c+p — |c_p 
is the asymmetry between the two minima, and /d -|- *7 = 2c!/c_e*^. 

We also calculate the statistical correlation function 

2 

Tfl 

{p{a 2 G 2 )p{aiGi))st = ^ [1 -hoi( 5 cosz/ti +/ 3 sinz/ti) 

-I-02 cos ij{t2 — ti)(cos vti — (3 sin uti) + 0102^ cos z/(t2 — H)]. ( 61 ) 

Again, we see that the fluctuations of the mass density are of the same order of magnitude 
as the expectation values. 


4. Probing a gravitational two-state (G2S) system 

Next, we present two different ways that mass density fluctuations of the two-level system 
described in Sec. 3 can be measured. In the hrst case, the probe, or the measuring apparatus, 
is a classical test-mass that records the Newtonian force continuously. In the second case, 
the probe is a quantum harmonic oscillator coupled to the gravitational two-state system 
described above. We consider the gravitationally induced effects to the oscillator’s evolution. 

4 . 1 . Classical test- mass probe: Fluctuating force 

Consider a test mass mo located near the conhning potential. We choose the coordinate 
system so that the two minima of the potential he at (±f, 0, 0) and the test mass probe at 
(0,*/,0), as in Fig. 1. 
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Figure 1. Force on a probe / detector exerted by a massive particle in a gravitational 
cat state between locations ±L/2 . 


We evaluate the x component of the Newtonian force exerted by the two-level system 
on the test particle. Again denote a = ±. If the system lies on the minimum of the potential 
at X = aL/2, the force Fx{a) exerted on the test particle in the x direction is 

GmmQL 


Fxia) = 


ro — —/ofl 


(62) 


2(9^ + LV4)3/2 

Assuming that the test mass is not allowed to move, the force Fx takes only two values /o 
and —/o- These values are correlated with the projectors Pa, Eq. (48). Thus F^ corresponds 
to a self-adjoint operator 

F = -foP+ + foP- = -/ods, (63) 


on the 2-state system’s Hilbert space. Thus, the gravitational force behaves as a quantum 
variable, its probabilities and correlations determined by quantum mechanics. 

Since Newton’s law of gravity is instantaneous, the test mass will be acted upon by 
the gravitational force at all times. Hence, the measurement output will be a time series of 
recorded force values. Thus, this set-up corresponds to a continuous measurement [16] of the 
G2S system. 

In Appendix A, we calculate the correlation functions of the force, obtained from the 
quantum probabilities for N successive measurements 


H(cii, ti, ci 2 ) ^ 2 j • ■ ■) (^Ntn) Tr FaNtM ■ • • ^a.2t2 PaitiPoPaitiPa2t2 ■ ■ ■ Pa„tn (64) 


Without loss of generality, we consider an initial |-t-) state. We hnd 


(E(t)) = 
{F(t')F(t)) = 


(65) 

( 66 ) 
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The decay constant T is defined as 


r = 


Z/V 


(67) 


where r is the temporal resolution of the probe. 

We see that both the mean value of the force and the correlations decay exponentially. 
However, the force fluctuations are highly non-Markovian, except for the limit u = 0, at 
which the force is trivially constant. The parameter T defines not only the decay rate, but 
also the memory time of the non-Markovian process. 


4-2. Quantum harmonic-oscillator probe: superpositions 

Now consider a quantum probe made of a harmonic oscillator of frequency u that is 
constrained to move along the horizontal axis as in Fig. 1. The Hamiltonian of the harmonic 
oscillator probe is 

Hp = ua^a (68) 

If the amplitude of the oscillations is much smaller than L, the length scale of the cat state, 
the force acted upon the oscillator along the x direction is approximately constant and equal 


to Eq. (62). This corresponds to an interaction Hamiltonian 

/o 


Hi = -foasx = 


:CT3(a + a^), 


(69) 


y/2moU 

Thus, the total Hamiltonian of the two state system interacting with the oscillator probe 


IS 


H = Hs + Hp Hj = u&i + ua^a + ga3{a + a^), (70) 

where the system Hamiltonian Hs is given by Eq. (|56| specialized to y = 0, i.e., it is 


equivalent to the Hamiltonian of a single-mode Jaynes-Cummings model with a coupling 
constant 

Often, the Jaynes-Gumming model is solved via the Rotating Wave Approximation 
(RWA). However, the RWA is a good only near resonance and for small values of the coupling. 
This is not the regime of physical relevance here. The harmonic oscillator is supposed to be 
a probe, a pointer variable in the sense of quantum measurement theory. For this reason, 
the coupling must be the dominant term in the total Hamiltonian, as seen, for example, in 
the standard model of a von Neumann measurement m- 

This implies that we must consider the deep strong coupling regime |1^ of the Jaynes- 
Cummings model that corresponds to g > u. While the model has been recently shown to 
be integrable and the spectrum of the Hamiltonian has been computed |19], it is still difficult 
to obtain closed expressions for the dynamics [50]. Here, we restrict to an exactly solvable 
regime which corresponds to vanishing tunneling rate and to the limit of weak tunneling 
rate which can be accessed perturbatively. This will serve as a first approximation for the 
problem under consideration. 
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4-2.1. Vanishing tunneling rate First, we assume that u « u, so that we can ignore the 
uai term in Eq. (56). We call the resultant Hamiltonian Hq: 

Ho = ujci^a + ga3{a + a^), (72) 

This approximation defines the adiabatic regime of the Jaynes-Cummings model, which 
is fully solvable. The evolution operator is 

^ A ( DHg/uj)e-^‘^^^^^Dig/u) 0 

where D{w) = stands for the displacement operator. 

As initial condition at t = 0, we assume that the oscillator is the vacuum state and the 
gravitational 2 state system is in a superposition state c+|-l-) J- c_|—). Evolving with (73), 
we find the state of the system at time t 

wICW) 


(73) 




c_ 


at)) 


where the path 


CW = --(l-e 

CO 


—iujt 


), 


describes an oscillation centered around Co = —- = 


fo 


corresponds to position xq = 


— fo 


0 2mou 


(74) 

(75) 

The center of the oscillation 




and momentum po = 0. 


The reduced density matrix for the oscillator degrees of freedom is 

pu^) = ic+nc(t))(c(t)i-7 ic_n - cw)(-cwi 
+ c\c2\at)){-at)\ + cic_\ - at)){at)\ 

The off-diagonal terms are distinguishable only if |((^o| — Co)P << 1 
g-4|Coh <-<- This implies that 

fo 


UJ" « 


mo 


mo- 


(76) 

or equivalently if 

(77) 


Eq. (77) is a necessary condition for the oscillator to act as a probe of the cat state. It 
demonstrates that the ultra-strong coupling limit of the Jaynes-Cummings model is the 
physically relevant regime. For sufficiently large mass mg, the oscillator may be a mesoscopic 
particle and its position at specific instant of time measurable. This would imply a resolution 
a for the position measurement smaller than the width xo- 

Hence, a position measurement at time t has probability equal to |c+p of finding the 
oscillator in the neighborhood of the phase space point ({t) and probability equal to |c_p of 
finding the oscillator in the neighborhood of the phase space point —({t). 

Of particular interest are the states |Co, +) = |Co) ® |+) and | — Co, —) = I — Co) ® |— ) 
where | ± Co) = -D(Co)|0) is a coherent state. We find that their time evolution involves only 
a time-dependent phase. 


—iHot I 


± Co. ±) ^ ± Co. ±), 


( 78 ) 


and thus they define stationary states of the system. 
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4-2.2. Rabi oscillations A finite value of u allows for transitions between the two 
gravitational quantum states, which induce transitions among the phase space paths of the 
oscillator. While the model is not exactly solvable, we can estimate the rate of such transitions 
using perturbation theory with respect to the tunneling rate u. In Appendix B, we show that 
to leading order in u, where 

cosnt — i sin i^t£)(2^o) 

—i sin utD{—2 (q) cos z/t 

As an estimate of the transition between the two gravitational quantum states, we 
compute the amplitude (—Co) +), between the stationary states |Co) +) and | —Co, —)• 

We hnd 

(-Co,-|Ot|Co,+) =-*sinz/t, (80) 

and thus the associated probability 

pit) = l(-Co,-|Ot|Co,+)P = sin^nt, (81) 

exhibits Rabi-type oscillations, with frequency u. 

5. Discussions 

The key findings of this work are enumerated in the Abstract. We make some additional 
comments below. 

1. Comparison with alternative quantum theories. At the Newtonian level, the gravitational 
field is completely slaved to the mass density. Poisson’s equation is the Newtonian limit of the 
Hamiltonian constraint of General Relativity and not a dynamical equation. Thus, once we 
have a theory that determines the dynamics and of mass density and its behavior under acts 
of measurement (in the standard quantum theory), we straightforwardly obtain a description 
for the Newtonian gravitational held. 

In particular, we need make no assumption that gravity is fundamentally quantum, 
because the quantization of gravity refers to the true degrees of freedom and not to the 
gauge-dependent variables that appear in the constraints. We only assume that there are 
no unknown physics in the relation between mass density and gravitational forces in the 
Newtonian regime. An example otherwise is the Newton-Schrodinger equation—see. Refs. 
[SHIS21ES] for the relevant discussion. 

Thus, the foundational question that can be settled by experiments invoking G2S systems 
is the following. Does the gravitational force remain slaved to the mass density even if 
the latter behaves quantum mechanically (i.e., when it fluctuates, or subjected to quantum 
measurements) ? 

The results derived here are based on standard quantum mechanics and (the 
nonrelativistic limit of) general relativity. Both theories have stood the test of time in their 
respective domains of validity. Results presented here could serve as a useful benchmark 
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for corresponding results from alternative quantum theories (AQT) to be compared and 
scrutinized. We expect discrepancies, even contradictions, appearing which should provide 
useful focal points for cross-examination of all AQTs in relation to the standard theories. 

2. Quantum Jumps of Spacetime Geometries. The cat state we consider here corresponds 
to the two minima of a potential, with a non-zero tunneling rate u between them. A 
transition between two orthogonal qubit states is standardly described as a quantum jump. 
The classic quantum jump experiments in atomic systems have shown that the duration of 
the jump is too small to be resolved [M], so jumps are effectively instantaneous. This implies 
the possibility of instantaneous jumps between two spacetime geometries that correspond 
to different mass distributions. The ’’instantaneous” part is not an issue here, because 
Newtonian gravity admits action at a distance, and in the non-relativistic regime we cannot 
explore issues of causality. However, the idea that quantum jumps can occur in the spacetime 
description because of the interaction of gravity with quantum matter is a new phenomenon 
of foundational value that is worthy of closer experimental exploration. 

3. From idealized models to realistic measurements. The two measurement schemes we 
have considered correspond to two major paradigms of quantum measurements, namely, a 
continuous measurement and a von Neumann measurement. Our results also exemplify 
typical behaviors of probes. We hud that a classical probe records a non-Markovian 
fluctuating force. A quantum probe interacting with the G2S system undergoes Rabi 
oscillations in the appropriate regime. For follow-up work along the present line of inquiry, 
state-of-the-art measurement schemes using atomic, molecular or light interferometric or 
electro- or opto-mechanical techniques would be very helpful. 
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Appendix A. Probabilities and correlations fnnctions for continuous 
measurement of the force 


We consider a continuous-time measurement of the operator (63) that corresponds to the 
gravitational force exerted upon a probe particle. 

Let T be the temporal resolution of the force measurements and let T stand for the 
duration of the measurement. We decompose the interval [0,T] into N = T/t subintervals, 
or equivalently A -|- 1 time steps, z = 0,1, 2,..., A. At each time step, there are two possible 
outcomes to the measurement, corresponding either to a = -|- or to a = — in Eq. (62). 
Without loss of generality we assume the initial state to be |-|-), i.e., c+ = 1, c_ = 0. 

A recorded series of values is represented by n-|-l integers ki that satisfy ~ ^+1, 

so that the hrst ko entries of the series are +, the following ki entries are —, and so on until 
the hnal kn entries that are (—1)”. The integer n < N stands for the number of jumps from 
-|- to — or the reverse. 

The probability associated with a sequence of A records is 

P{ai, ti] 02 ,^ 2 ) • • •; = Tr -Pajvtjv • • • ^a2t2^aitiP0-Paiti-Pa2t2 ■ ■ ■ ^a„tn ■> (-^T) 

where a* = ±. 

Taking times ti = ir, i = 1,... N and representing the sequence of records by the n 
integers fcj, we obtain 

F(fcl, ^2, . . . , kn) = K + |f>.|+)P^'^-")|(-|f>.|+)r", (A.2) 

where we have used the fact that (+|Fr|+) = {~\Ut\—) in a symmetric conhning potential. 


Remarkably, the probability (A.2) depends only on the number n of jumps. 


Since |(+|Lh|-|-)f = cos^ ^ and |(—|Lh|-|-)f = sim Eq. (A.2) becomes 


/ i/r\2(A^-n) 

P{ki, ^ 2 , • • • , kn) = [cos — 1 




Z/r\2ri 

smyj 


In the regime z/r << 1, 

F(fcl,/C2,- 

where 


.,kn) = X^e-^^^-^\ 


A = 


Z/V2 


(A.3) 


(A.4) 


(A.5) 


In order to hnd the correlation functions associated with the probabilities (A.4), we compute 


the conditional probability F(a 2 , Z 2 |aiAi) that for 02 at time step Z 2 given a measurement 


outcome ai at time step R. This equals the sum of the probabilities (A.4) over all sequences 


Oi over time steps z, such that A < z < Z 2 , such that = oi and 0*2 = 02 


To calculate this we hrst note that the probabilities (A.4) do not depend on any initial 


moment of time, so the conditional probability is a function of m = A ~ A- Thus, we need 
to compute a function g{a 2 ,ai,m) = F(a 2 , Alcq, A)- 
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First, consider oi = 02 = +. Only seqnences with an even nnmber n = 2/ of jnmps 
contribnte. The nnmber of different seqnences of length m + 1 that start from + and contain 
n jnmps is (™). Hence, 


72] 


9(+,+;m)=(cos^)^”5; 


1=0 


- sm z/r 
2 


4 / 


Similarly, we find 


= (cos y) 


2771 


[(m-l)/ 2 ] 

E 

1=0 


m 

21 + 1 


1 

- sin VT 
2 


4Z+2 


We employ the binomial identity (1 + x)™' = (7)^^ obtain 

^(+,+;m) = ^ (cos 



1 1 

('“y) 

(1 -7 - sin^ z/r)”" -h (1 - - sin^ z/r)"* 

/ z/r\2»". 


(““y) 

1 + - sin^ z/r ”* - 1 - - sin^ z/r W 
'4 4 


(A.6) 

(A. 7) 

(A.8) 

(A.9) 


By symmetry, g{-, m) = g{+, +; m) and g{+, m) = g{-, +; m). 

We assnme an initial state |+). The expectation valne of the force at the time-step m is 

1 


{F{m)) = -fo[g{+, +; m) - g{-, +; m)] = -/o 
The correlation fnnction at time steps mi and m 2 > mi is 


2 X 2 

cos y(l —-sin z/r) 


(A.IO) 


{F(m 2 )F(mi)) = [g{+, -F; m 2 - mi) - g{-, -F; m 2 - mi)] 

X [g{+, -h; mi) -h g{-, -F; mi)] 

= /o (cos^ YJ “ 4 (1 + - sin^ z/r)”"i (ATI) 

In the regime z/r << 1, 


(F(m)) = -/oe-2"™ 
{F(m2)F(mi)) = ^ 


where A is given by Eq. (A.5). 

The continnons-time limit follows by setting m 
T = 2A/r = Then, Eqs. (65, 66) follow. 


(A.12) 

(A.13) 


tjr and dehning the decay constant 












Probing a Gravitational Cat State 

Appendix B. Perturbative propagator for the Jaynes-Cumming model 


23 


We compute the evolution operator e where H is given by Eq. (70). We write 
H = Ho + V, where Ho is given by Eq. (72) and V = uai. We work in the interaction 

(B.l) 


representation. We dehne the unitary operator 


— pHot -iHt 


which is a solution to the equation 

'wfi‘ = 


with initial condition Oo = 1. In Eq. (B.2) 
V{t) = = 

in terms of 


0 St 

S* 0 


St = £>(C„)e"»‘C(-2Co)e-"»‘£>(Co) 

where D{z) is the standard displacement operator. 

To leading order in perturbation theory, the solution Ot can be expressed as 

O, = Jo' 


(B.2) 


(B.3) 


(B.4) 


(B.5) 


i.e., the time-ordered exponential that solves Eq. (B.2) is substituted by an ordinary 
exponential. We hnd 


dsS{s) = D{(a) £)(Co). 


(B.6) 


For ut » 1, the time integration over a periodic function is negligible, so 

/g ~ t. Hence, J^dsS{s) = i){2(o), and 


dsV{s) = vt 


0 D{2Co) 
D{-2Co) 0 


(B.7) 


By Eq. (B.5), we obtain Eq. (79) for Ot- 
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